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ABSTRACT 

X-ray observations of hot gas in galaxy clusters often show steeper temperature gradi- 
ents across cold fronts - contact discontinuities, driven by the differential gas motions. 
These sharp (few kpc wide) surface brightness/temperature discontinuities would be 
quickly smeared out by the electron thermal conduction in unmagnetized plasma, sug- 
gesting significant suppression of the heat flow across the discontinuities. In fact, the 
character of the gas flow near cold fronts is favorable for suppression of conduction 
by aligning magnetic field lines along the discontinuities. We argue that a similar 
mechanism is operating in the bulk of the gas. Generic 3D random isotropic and in- 
compressible motions increase the temperature gradients (in some places) and at the 
same time suppress the conduction by aligning the magnetic field lines perpendicular 
to the temperature gradient. We show that the suppression of the effective conduc- 
tivity in the bulk of the gas can be linked to the increase of the frozen magnetic field 
energy density. On average the rate of decay of the temperature fluctuations d{ST^)/dt 
decreases as (_B^)^^/^. 
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1 INTRODUCTION 

X-ray observations of galaxy clusters reveal significant spa- 
tial fluctuations of the gas temperature in a range of spatial 
scales (e.g., Markevitch et al. 2003). Given a temperature 
map with prominent fluctuations, it is possible to calculate 
an upper limit on the effective thermal conductivity, pro- 
vided that the lifetime of the fluctuations can be estimated. 
It turns out to be at least an order of magnitude lower than 
the Spitzer conductivity for unmagnetized plasma (Ettori & 
Fabian 2000; Markevitch et al. 2003). 

Heat conduction in the intracluster medium (ICM) is 
primarily along the field lines because the Larmor radius of 
the particles is very small compared to the collisional mean 
free path (Braginskii 1965). The ICM undergoes turbulent 
motion in a range of spatial scales (Inogamov & Sunyaev 
2003; Schuecker et al. 2004; Schekochihin & Cowley 2006; 
Subramanian et al. 2006; Zhuravleva et al. 2011). As the 
magnetic field is, to a good approximation, frozen into the 
ICM, the fleld lines become tangled by gas motions and their 
topology changes constantly. To study thermal conduction 
in such tangled magnetic flelds, some authors employed the 
random- walk representation of the magnetic fleld lines (Trib- 
ble 1989; Malyshkin 2001). Four main effects should be con- 



sidered. Firstly, parallel thermal conduction along stochas- 
tic magnetic field lines may be reduced because the heat- 
conducting electrons become trapped and detrapped be- 
tween regions of strong magnetic fleld (magnetic mirrors, see 
Chandran & Cowley 1998; Chandran et al. 1999; Malyshkin 
& Kulsrud 2001; Albright et al. 2001). Secondly, diffusion 
in the transverse direction may be boosted due to spatial 
divergence of the fleld lines (Skilling et al. 1974; Rechester 
& Rosenbluth 1978; Chandran & Cowley 1998; Narayan & 
Medvedev 2001; Chandran & Maron 2004). Thirdly, there is 
effective diffusion due to temporal change in the magnetic 
field ( "fleld-line wandering" ) . Finally, if one is interested in 
temperature fluctuations and their diffusion, one must be 
mindful of the fact that the temporal evolution of the mag- 
netic field is correlated with the evolution of the temper- 
ature field because the fleld lines and the temperature are 
advected by the same turbulent velocity field. 

In this paper, we focus on the last effect. The more con- 
ventional approach, often used to estimate the relaxation of 
the temperature gradients, is to consider the temperature 
distribution as given and study the effect of a tangled mag- 
netic fleld on the heat conduction. However, the direction 
and value of the fluctuating temperature gradients are not 
statistically independent of the direction of the magnetic 
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field lines because the latter are also correlated with the tur- 
bulent motions of the medium. We argue that, dynamically, 
the fluctuating gradients tend to be oriented perpendicular 
to the field lines and so heat fluxes are the more heavily 
suppressed the stronger the thermal gradients are. We also 
establish the relationship between the average conductivity 
and the growth of the magnetic energy density. 

The structure of the paper is as follows. In § 2, we pro- 
vide a qualitative explanation of the correlation between the 
temperature gradients and the magnetic field direction, ac- 
companied by a number of numerical examples. In § 3, a 
theoretical framework for modelling this effect is presented 
and the joint PDF of the thermal gradients, the angles be- 
tween these gradients and the magnetic field lines and the 
magnetic-field strength is derived in the solvable case of a 
simple model velocity field. The connection between the ef- 
fective conductivity and the increase of the magnetic energy 
density is established. Analytical results are supplemented 
by numerical calculations in § 3.4, which extrapolate our re- 
sults to the case of a more general velocity field. In § 4, we 
provide a general discussion (including of the limitations of 
our treatment). Finally, in § 5, we sum up our findings. 



2 QUALITATIVE DISCUSSION 

We consider a volume of plasma with high electric conduc- 
tivity and frozen-in magnetic field tangled on a scale much 
greater than the mean free path of the particles. We also 
assume the plasma motions to be incompressible, which is a 
good approximation for subsonic dynamics. Across the pa- 
per we treat the temperature as a passive scalar. 



2.1 Illustrative example: conduction bet^veen 
converging layers of magnetised plasma 

Consider two parallel layers of an incompressible medium 
vertically separated by distance h with temperatures Ti 7^ 
T2. This is illustrated in Fig. 1: the direction of the field line 
is shown with the inclined solid line, making an angle 6 with 
the vertical, so cos6' — h/\'hP + P, where I is the horizontal 
distance between the footpoints of the field line anchored 
in the two layers. An incompressible flow with dyUy < 
reduces h and increases I so that / x ^ is conserved (in the 
absence of tangential shear). Here we are interested in the 
heat exchange between the layers, i.e., only the component 
of the heat flux along the temperature gradient Qvt has to 
be calculated; 



OvT = x(b • VT) cos 6* = X 



T2-T1 



h 



VWTP \/h?TW 



Let h{t) = /io/(i) and l{t) = k)/f{t). Then 



!VT = X- 



T2-T1 



f 



ho 



f + {lo/hoff-^' 



(1) 



(2) 



where x is the parallel thermal diffusivity coefficient (Bra- 
ginskii 1965), which is assumed constant across the volume 
for simplicity. Therefore, in the limit of / — > 0, Qvt — >■ if 
lo 7^ 0. Similarly, when / — > cx), Qvt — > 0. The decrease of 
the heat fiux at / > 1 is simply due to the increase of the 
distance between the plates and corresponding decrease of 
the temperature gradient. The decrease at / < 1 is due to 




Figure 1. Correlated changes of the temperature gradients and 
the inclination of the magnetic field lines in the ease of a converg- 
ing incompressible flow: plane parallel layers at different temper- 
atures. Converging flow with dyUy < reduces h and increases 
the temperature gradient (T2 — T\)/h, but suppresses heat flux 
(Fig. 2). The solid line represents the direction of the magnetic 
field. If the medium is incompressible then I X h is conserved (in 
the absence of tangential shears) . 




Figure 2. Suppression of the heat flux along the temperature 
gradient between two approaching/receding plates as a function 
of distance / between the plates, when the medium between the 
plates is threaded by tangled magnetic field [see equation (3)]. 
At the initial moment (/ = 1), all angles between the magnetic 
field direction and the plates are equally probable. The decrease 
of the heat flux at / > 1 is simply due to the increase of the 
distance between the plates and corresponding decrease of the 
temperature gradient. The decrease at / < 1 is due to systematic 
increase of the angle between the field lines and the direction of 
the temperature gradient. 
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systematic increase of the angle between the field lines and 
the direction of the temperature gradient. 

If at some moment the field lines are tangled in such a 
way that all angles 6 are equally probable, then parametriz- 
ing compression/stretching along y by the same factor / and 
averaging over 6 gives us the suppressed heat flux along the 
temperature gradient: 



>VT = X- 



Ta 



2/ 



fto 



(see Fig. 2). 



(3) 



/2 + 1 

Thus, increasing the temperature gradient by squeezing the 
layers of the gas does not boost the heat exchange between 
them but rather makes it smaller. A qualitatively similar 
situation might occur at the cold fronts - contact discon- 
tinuities formed by differential gas motions, a very simple 
model of which is discussed in the next subsection. 



2.2 Astrophysical example: model of a cold front 

Chandra observations of galaxy clusters often show sharp 
discontinuities in the surface brightness of the ICM emis- 
sion (see review by Markevitch & Vikhlinin 2007). Most of 
these structures have lower-temperature gas on the brighter 
(higher-density) side of the discontinuity, suggesting that 
they are contact discontinuities rather than shocks. In the 
literature, these structures are called "cold fronts" . Because 
of the sharp temperature gradients, the limits on the thermal 
conduction derived for the observed cold fronts are strong 
(see, e.g., Ettori & Fabian 2000; Vikhlinin et al. 2001; Xiang 
et al. 2007). 

In the majority of theoretical models, the formation of 
a cold front involves relative motion of cold and hot gases. 
Here we consider the case of a hot gas flowing around a 
colder, gravitationally bound gas cloud, which is a proto- 
typical model of a cold front. For simplicity, we assume that 
the velocity field can be approximated with a 2D potential 
flow past a cylinder, while the initial temperature is sym- 
metric around the cylinder. The initial temperature distri- 
bution and stream lines of the flow are shown in the left 
panel of Fig. 3. The middle panel shows the field lines of a 
random magnetic field superimposed on the initial temper- 
ature distribution. The evolved temperature and magnetic 
field are shown in the right panel of Fig. 3. Stretching of 
the fluid elements near the stagnation point along the front 
leads to the contraction of the same elements in the direc- 
tion perpendicular to the front. This configuration has been 
considered in a number of studies of the cold fronts (see, e.g., 
Asai et al. 2007; Churazov & Inogamov 2004; Roediger et al. 
2011). Qualitatively, it corresponds to the situation sketched 
in § 2.1 and Fig. 1, which naturally leads to the field lines 
orthogonal to the temperature gradient at the front. 




Figure 3. Alignment of the field lines perpendicular to the tem- 
perature gradient for the velocity field characteristic of a cold 
front. A potential flow past a cylinder is used in this example. 
The left panel shows the initial temperature distribution (color 
image) and stream lines of the velocity field. The middle panel 
shows a random tangled magnetic field superposed on this initial 
situation. The right panel shows the time-evolved temperature 
map and magnetic field lines (superposed contours) in such a flow. 
The flow boosts the temperature gradient at the cold front and 
at the same time stretches the field lines along the lines of con- 
stant temperature. In the resulting configuration, the field lines 
are essentially perpendicular to the sharp temperature gradient 
at the front. 



dT 
dt 



0, 



(5) 



where B is the magnetic field, u the velocity field, T tem- 
perature and d/di = d/dt+u-V . We have neglected thermal 
and magnetic diffusivities. Let g be the unit vector in the di- 
rection of the temperature gradient, b the unit vector in the 
direction of the field line, B the magnetic field magnitude 
and G the temperature gradient magnitude, so B = Bb, 
Vr = Gg. The above equations imply: 

-g = -Gg-(Vn).g, (6) 

AB 

m 

d/i 

dt 

where /j. — b ■ g, the cosine of the angle between B and 
vr. From these equations, we can immediately infer the 
following equation for b ■ \7T = iiG, a quantity proportional 
to the parallel heat fiux: 

dln(^G) _ dhiB 

At ^ At ■ ^> 

Thus, locally, the heat flux decreases as the fleld strength 
grows. 



= Bb-{Vu)-b, 

= fi[g-(\/u)-g-b-{\/u)-b], 



(7) 
(8) 



2.3 Local correlation between the magnetic field 
strength and the heat flux 

Let us now discuss the suppression of the local heat flux in 
more general terms. Consider the induction equation for an 
incompressible medium and the advection equation for the 
temperature: 

AB 



= BVu, 



(4) 



2.4 Numerical example: a random 2D velocity 
field 

In this example, we consider a random temperature distribu- 
tion and a random magnetic fleld in a random (5-correlated 
Gaussian incompressible 2D velocity field (Fig. 4). The tem- 
perature T(x, y), the magnetic field B{x, y) and the velocity 
field u{x, y) (assumed incompressible, V • u = 0) are mod- 
eled as superpositions of Fourier harmonics with random 
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phases and amplitudes. The temperature and the magnetic 
field are advected according to equations (4) and (5). The 
velocity field is varied at each time step as a (5-correlated 
in time field. The initial conditions are shown in the top 
panel of Fig. 4; there is no initial correlation between the 
temperature gradients and the orientation of the field lines. 
With time, preferential stretching/squeezing of the fluid el- 
ements leads to alignment of the field lines along the iso- 
temperature lines (see bottom panel in Fig. 4). This happens 
in all regions where the stretching/squeezing is sufficiently 
strong. As a result, the field lines are mostly perpendicular 
to the direction of the temperature gradient in all regions 
where the gradient is large. Intuitively, one expects this ten- 
dency of local alignment between the magnetic filed and the 
isotherms to manifest itself statistically, in a turbulent con- 
ducting medium. In the next section, we work out a simple 
statistical model of this process. 



3 HEAT CONDUCTION IN A STOCHASTIC 
VELOCITY FIELD 

Here we treat the suppression of the heat conduction using 
an analytically solvable model that allows us to predict the 
statistical distribution of the cosine of the angle between the 
thermal gradient and the field line (fi) , the magnitude of the 
thermal gradient (G) and the magnetic-field strength (B). 
After the joint probability distribution function (PDF) of /i, 
G and B is derived (§ 3.5), we will be in a position to assess 
how statistically prevalent the behaviour discussed in § 2.4 
is, but we will preface this detailed calculation with some 
simpler arguments to quantify the suppression of the heat 
flux. 




3.1 Relaxation of temperature fluctuations 

Let us restore heat conduction in equation (5): 

^=V-(xbb-Vr), (10) 

where x is the parallel thermal diffusivity coefficient (Bra- 
ginskii 1965). Then the volume- averaged rate of change of 
the rms temperature fluctuations is 

d(<5r^} 



df 



= -2x(|6-V5r|=') = -2x(M'G'). 



(11) 



Thus, the average value of fi^G^ characterizes the rate at 
which local temperature variations are wiped out by the 
thermal conduction. 



3.2 Kazantsev-Kraichnan model 

We consider the magnetic field to be so weak that it does 
not affect the velocity field. This condition is only satisfied if 
the magnetic energy density is much lower than the kinetic 
energy density of the plasma motions. This means that our 
model does not describe the saturated state, when these 
energy densities become comparable. The non-saturated 
regime could be a common transient situation in the ICM, 
at least locally, in the sense that at any given time, the mag- 
netic field is amplified up to the saturation value only in a 
small fraction of the volume. 



Figure 4. Alignment of the field lines perpendicular to the tem- 
perature gradient for a stochastic (5-correlated Gaussian incom- 
pressible velocity field, modeled as a superposition of Fourier har- 
monics with random phases and amplitudes. The top panel shows 
the initial random temperature distribution (color) with the field 
lines of a random magnetic field superposed (they are uncorre- 
lated with temperature). The bottom panel shows the same fields 
later on in the evolution. In the evolved image, field lines fol- 
low the lines of constant temperature in the regions where the 
temperature gradient is large. 



We will wish to calculate the joint PDF p{fj,,G,B;t), 
where /i and G are defined in § 2, and investigate the evo- 
lution of the relevant correlations, viz., (jj^G^) (see § 3.1). 
To do that, we need to average the dynamical equations for 
g, b, G and B over all realizations of the stochastic velocity 
field. The equations are: 



di 

H 
dG 
dt 
dS 

di 



-{Om - 9 9 )g OmU , 



/ ck J k ] i\j 7n r^ i 

(di — b b )b dmU , 



— — = —Gg g dmU , 



(12) 
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where summation over repeated indices is implied. 

This problem is solvable analytically for a Gaussian 
white-in-time velocity field (Kazantsev f968): 



{u*(f, x)u' {t ,x )) = 5{t — t )n-' {x — x), 



(13) 



where k'-' is the correlation tensor, whose form can be deter- 
mined from symmetry and incompressibility considerations. 
We may assume the medium to be isotropic and homoge- 
neous. Consider a spatial scale much smaller than that of 
the stochastic velocity field. Then, at any arbitrary point, 
the velocity can be expanded in linear approximation: 



u*(f, x) — al^{t)x" 



(14) 



where (T^(i) — dmu" and we have assumed M*(t, 0) = with- 
out loss of generality (otherwise change the reference frame). 
Then the velocity gradients satisfy 



(^(*'")^('''"')) 



{aUt)<Ti{t')) 



5{t - t')K, 



mn 1 



(15) 



where 



d't^'^y) 



dymdyn 






J/=0 



where k = 1/ Teddy, Teddy is the turn-over time of the turbu- 
lent eddies and 

1 



5''5„ 



d+1 






(16) 



is the inevitable tensor form of K^^n for an isotropic incom- 
pressible medium of dimension d (= 2, 3). This is the so- 
called Kazantsev-Kraichnan model, which has been a pop- 
ular tool for modelling the properties of small-scale dy- 
namo and passive-scalar advection in turbulent media (e.g., 
Chertkov et al. 1999; Balkovsky & Fouxon 1999; Boldyrev & 
Schekochihin 2000; Schekochihin et al. 2002, 2004; Boldyrev 
& Cattaneo 2004, and references therein). 



3.3 Relation between magnetic field amplification 
and suppression of conduction for the 
white-in-time velocity field 

Before presenting the full statistical calculation, we wish to 
give a relatively simple one that establishes the connection 
between the relaxation rate of the temperature fluctuations 
and the magnetic energy density. The heat flux along the 
held line /iG is inversly proportional to the length of a fleld- 
line segment s. Therefore, one can relate the change of the 
mean square heat flux {y^G^), which is also the decay rate 
of the temperature fluctuations (see §3.1), to the growth of 
the magnetic-energy density as follows: 



{B^) cc (s^), i^^G^) cc {1/s^). 



(17) 



As explained in § 3.2, we assume an isotropic linear 
random velocity held. Let it be piecewise constant in time 
over intervals Tc and completely uncorrelated for Ai > Tc- 
Assume further that the amount of stretching of any fluid 
element over individual time intervals of duration ~ Tc is 
small compared to the size of the element, which amounts 
to a model of white-noise held. Under these assumptions. 



it is easy to obtain the PDF of s as a function of time t 
in the limit I/tc 3> 1. The evolution of each component of 
the separation vector x of any two locations frozen into a 
velocity field constant over time interval Tc is 



1 



(18) 



x\Te) ~ x\Q) + Teo]x'{Q) + lr»^a;"(0) + 0{T'e), 

where o^j is the velocity gradients matrix [see equation (14)]. 
Since we are dealing with a random isotropic field, we can 
set a3(0) = (1, 0, 0) at i = 0. Then 



X {Tc) ^ 1 + TcO-i + -Te a^a^ 
2; (Tc) = Tca{ + Oij'e). 



o{tI), 



(19) 



We are interested in the time evolution of the "stretching 
factor" s^ — \x\^ . For one "act of stretching" , equation (19) 
implies 

Inspire) = 2Teo\-2T'^e{(^l? + rl(j{a{+Tla]a{+0{Tl).{2Q) 

For t ^ Tc, the calculation of s^(t) reduces to summation of 
TV = t/Tc ^ 1 such independent stretching episodes: 



In s (t) — y^ In s {tc 



(21) 



After applying the central limit theorem to 
(1/iV) J2 In s^(t-c), one readily gets the PDF of s^: 



Pi^') = ^ 



1 



_(,„. 



^r 



V2^ 



(22) 



where 



= 2. Til 



'^eddy 
d 

2Tii + J2 (n 



rjiil 
' J-li 



t 



'^eddy 



(23) 



where Teddy and T^„ are defined at the end of § 3.2. We have 
taken 5(0) = 1/tc in equation (15). Using equation (17), we 
get 



(B^) oc e 



"ls+CTj/2 



{^^'G^} 



(24) 



This leads to a simple relation between the growing 
magnetic-energy density and the evolution of the mean 
square heat flux: 



{^ G ) oc {B Y , whore p 



-nis + as/2 
ms + al/2 ■ 



(25) 



For an incompressible velocity field in 3D, using equa- 
tion (16), we get p — —1/5. This is a statistical version 
of the dynamical equation (9). It implies that on average, as 
the magnetic-energy density grows, the rate of decay of the 
temperature fiuctuations is reduced, although the efficiency 
of this reduction is modest (p is low). This is because {^G^) 
is dominated by regions of low stretching while (B^) by re- 
gions of high stretching [equation (17)] and the distribution 
of these is highly intermittent. 



3.4 Finite-time correlated velocity field 

How sensitive is this result to the obviously unphysical as- 
sumption of zero correlation time? Here, we numerically cal- 
culate the PDF of s in a random incompressible 3D velocity 
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field evolving according to a Langevin equation with a finite 
correlation time. This is a generalization of the (5-correlated 
case considered in § 3.3. 

We consider a large number of independent field-line 
segments, each one placed in its own stochastic incompress- 
ible velocity field, given by equation (14), where the velocity 
gradient satisfies 



dt 



-a^ + 9mffl\ 



(26) 



where Tc is the correlation time and the gradient of the 
stochastic acceleration o* satisfies 



{9™a'(t)9„a^(t')) = <5(i - tV'r, 



mn • 



(27) 



Here A^ is the noise amplitude and the dimensionless tensor 
T^„ is fixed by isotropy and incompressibility as given by 
equation (16). It is possible to define the efi'ective turn-over 
time of turbulent eddies Teddy in much the same way as we 
we did for the (5-correlated case: 



{aUO)<ri{t))dt^^A\!T};i„ 



1 



2r< 



eddy 



(28) 



where the exact solution of the Langevin equation (26) has 
been substituted. Thus, r^ddy ~ 1/{tcA)'^. 

In view of equation (17), the evolution of {y?G^) and 
{B^) can be easily calculated from the distribution of the 
segment lengths. Here we do this for a range of values of 
the ratio rj — Tc/T^ddy In § 3.3, we treated the case rj — > 
analitically, whereas for a physically sound case, 77 ~ 1 
because typically turbulent velocities decorrelate over their 
eddy turn-over times and fiuid elements are stretched by 
order-unity amounts over the same time scales. The results 
are shown in Fig. 5. Even though the growth/decay rates 
of (-B^) and (/I'^G^) do change with correlation time, their 
relative behaviour appears to be invariant, viz.. 



{^i^G'^) ex {bY^^\ 



(29) 



practically the same as for the (5-correlated regime [cf. equa- 
tion (25)] 

Thus, finite correlation times do not change the form 
of the effective conduction-magnetic-energy-density relation, 
only modifying the time dependence. This result gives us 
some confidence in the Kazantsev-Kraichan velocity as a 
credible modelling choice. 



3.5 Statistics of the heat flux 

In this section we will finally derive the full joint statistical 
distribution of the magnetic field and the temperature gra- 
dient and hence the detailed correlations between the heat 
flux, the fleld strength and the relative direction of the mag- 
netic fleld and the temperature gradient. 

For a velocity fleld given by equation (14), we can write 
equations (12) for g, b, G and B as follows: 

tg = -yo^ - g g )g <^m, 

d] k / ck J kj i\ ] m i 

tb = [Si -b b )b cr^, 

dtG = -Gg^g'^a]^, 

dtB = Bb'b'^al,. (30) 

There are no advection terms here due to the homogeneity 
of the gas (so we can consider equation (12) at cc = 0). 



The details of the derivation of the equation for the joint 
PDF p{iJ.,G,B;t) are presented in Appendix A. The result 

is 



dtp ■ 



2{d- 



2rf(l - fi ){)idfi)idfi - daGfid^ - dsBfid^) 



+{d - iKdoGdoG + ObBObB) + 2(1 - ti'd)dGGdBB 
+d{d + 1 - 2dn'^){2fid^, - dcG - dsB) 



+2d\l-dfi-' 



P, 



(31) 



where d is the dimension of space. From now on, we only 
consider d = 3. 

Multiplying both sides of equation (31) by /i^G^ and 
integrating, we find 



^t{^^^G^) = -^(m'g^), 



(32) 



so the mean square heat flux decays exponentially in time. 
Then, recalling equation (11) for the rate of smoothing of 
the temperature fluctuations. 



-Kt/2 



0. 



(33) 



d(ST^) 
-^ — oc 
dt 

We observe that the relaxation rate of the temperature fluc- 
tuations decreases signiflcantly on time scales of the order 
of the turn-over time of the turbulent eddies {k = l/T^ddy)- 
It is also possible to reproduce the relation for the mean 
square heat flux as a function of the magnetic energy den- 
sity [equation (25)]. Multiplying equation (31) by B^ and 
integrating, we obtain the evolution of the magnetic energy 
density: 

dt{B'} = ^n{B'}. (34) 

This result, combined with equation (32), leads to the rela- 
tion established in 5 3.3: 



(m'g^) = {bY'^' 



(35) 



We expect that the temperature gradients and the mag- 
netic fleld lines will become perpendicular to each other. Let 
us then flrst investigate the limit of /i -> 0, in which equa- 
tion (31) can be solved analytically. Let a; = In/i, y — InG 
and z = In B. Then the joint PDF of these variables is 
h{x,y,z;t) = p{fi{x),G{y),B{z);t)e'^'^^^^, where the last 
factor is the Jacobian of the transformation of variables. 
Taking /i -> in equation (31), we find that h satisfies 



dth 



'j'^xx \ '^yy I i^zz ^\'^xy ~T~ •^xz) \ ily 



+Z{2hx — hy — hz 

Let us now write h in the following form: 

h{x, y, z; t) = f(x, y; t)S{x + y + z). 



(36) 



(37) 



Substituting this ansatz into equation (36) , we find that the 
factorization goes through and / satisfies 



dtf 



•^Jxx + Jyy '^Jxy + o(2/x Jy) 



(38) 
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Figure 5. The decrease of the mean square heat flux {fi'^G'^) for the time-correlated velocity field and different ratios 77 = Tc/r^ddy 
(numerical results). While growth/decay rates of (B^) and [fj^G"^) change with correlation time, their relative behaviour is practically 
the same: (m^G^) oc (B2)-0.2^ 



This factorization implies that in the limit /i — > 0, /iG oc 1/ B 
independently from the initial conditions. This result was 
anticipated in § 3.3, where we took the ratio of pG and 1/B 
to be the same for all the segments of the field lines at the 
initial moment. 

Let us make another transformation: ^ = a; = In/i 
and r] = X + 2y — ln(/iG^) to separate variables in equa- 
tion (38). The joint PDF of these two variables, w{£_, -q; t) = 
f[x[£,),y{£,,ini);t), satisfies 



dtw = - ( 3it)^5 -I- -u;^,, 



6«;j 



(39) 



This equation can be readily solved 
1 



w{i,T];t) = 



VoTTKf 

rt('7-'7')' 



de'dr,7(e',r,';0)e-^[*'''+<«-«''] 

(40) 



Notice that along with diffusion in both variables, the PDF 
drifts to ^ — > —00, i.e, to smaller /i. So there is a continued 
tendency towards mutually perpendicular orientation of the 
thermal gradients and the field lines. 

If one is interested how the joint PDF of /i and G be- 
haves in the case of ^i order unity, full equation (31) inte- 
grated over the magnetic field strength has to be solved. 
Technically speaking, we are obliged to do this in order 
to ascertain that the limit ^ — >■ was the relevant one 
to consider, i.e., that the joint distribution of ^l and G 
moves towards smaller ^ independently of initial condi- 
tions. Again, to separate variables, we employ the variables 
^ = In/i and 77 = ln(/iG^). The PDF of these variables. 



w(C,');i) = /p(M(C),G(C,M),-B;i)e^«+'')dB, satisfies 



dtw = 



3(l-e''*)to^5 + (l-f3e^'=)m 



2«>, 



+6(1-26^^)1^5 -12e-"'w 



2«, 



(41) 



In order to solve this equation numerically, it is convenient 
to rewrite it in the divergence form as follows: 

a,«;= j|95[2(l-e^«) + (l-e^«)ad 



+dr,{l + ie^)dr,\w 



(42) 



Numerical solution of this equation is presented in Fig. 6. 
With time, the maximum of the PDF does indeed shift to- 
wards smaller ^, demonstrating that the temperature gra- 
dient and the magnetic field vector are becoming ever more 
orthogonal to each other. One can replot this graph in co- 
ordinates /iG (heat flux) and G to observe that the rate of 
smearing of the temperature fluctuations in equation (11) is 
correlated with the magnitude of the temperature gradients 
(Fig. 7) in such a way that sharper gradients on average 
tend to be wiped out slower due to smaller corresponding 
values of /iG. 



4 DISCUSSION 

We now briefly discuss the assumptions made in our model 
and its implications. 

1) The ordering of scales in the problem considered in 
this paper obeys the following relations: 



© 0000 RAS, MNRAS 000, 000-000 



S. V. Komarov, E. M. Churazov and A. A. Schekochihin 




Figure 6. Evolution of the joint PDF of fj, and G at regular time intervals from t = to t = T^ddy (turn-over time of the turbulent eddies) 
obtained via numerical solution of the equation (42). The maximum of the function drifts to the region where the thermal gradients and 
the field lines are perpendicular (^ — > 0). 
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Figure 7. Evolution of the joint PDF in terms of heat flux /iG = \b- VT| and G = |Vr| at the same times as in Fig. 6. Sharper gradients 
tend to be wiped out slower due to the smaller corresponding values of the heat flux. 



pe ^ ^mfp ^ ^ ^ -^M 



(43) 



where I is the characteristic size of the region we deal with, 
Pe is the electron Larmor radius, Xmfp is the electron mean 
free path and At, is the typical size of a turbulent eddy. 

The limit I <C A„ simpUfies the calculation of the field- 
line stretching because the linear expansion of the velocity 
field can be used [equation (14)]. This allows for analytic 
treatment of the problem. Note that kinematic dynamo nat- 
urally sets the parallel correlation length of the magnetic 
field Aflii to be ~ A„ (Schekochihin et al. 2002, 2004). 

The condition Amp/ ^ I allows us to apply the thermal 
conduction equation (10) at these spatial scales. Due to the 
fact that in the kinematic-dynamo regime, Au ~ -^s||i ^^ 
also have Xmfp ^ ^b\\- This limit being assumed, we can 
ignore the magnetic mirroring effects because the electrons 
are free to escape magnetic traps by changing their pitch 
angles by collisions (Chandran & Cowley 1998; Chandran 
et al. 1999). 

A typical value of Xmfp ~ 

20kpc(T/10* K)^(ne/10"^ cm"^)"^ due to Coulomb 
collisions varies in the cluster cores from 0.01 to 20 kpc 
depending on temperature and density. For example, in 
the core of the Coma cluster the mean free path is ~ 5 
kpc (Churazov et al. 2012); in M87/Virgo it is much 
smaller, viz., Xmfp ~ 0.01 kpc, due to lower temperature 
and higher density (Churazov et al. 2008). On the other 



hand, the value of At, can be in the range of 10 kpc to 
200 kpc (Inogamov & Sunyaev 2003; Schuecker et al. 2004; 
Schekochihin & Cowley 2006; Subramanian et al. 2006; 
Zhuravleva et al. 2011; Kunz et al. 2011). Therefore, our 
analysis is relevant for scales between ~ 0.1 to ~ 100 kpc. 
Some of these scales are directly resolvable with Chandra or 
XMM-Newton, suggesting that the observed substructures 
in the temperature maps should have the magnetic field 
lines roughly aligned with the iso-temperature contours. 

2) The assumption of incompressibility (needed to use 
equation (16) for the description of the velocity field) is valid 
as long as the gas velocities are subsonic. This is reasonable 
for the ICM, except for cases of strong mergers or AGN- 
driven strong shocks in the very core of a cluster. The com- 
parison of cluster mass estimates from X-ray data and lens- 
ing or stellar kinematics (e.g., Churazov et al. 2008) and 
simulations (e.g., Lau et al. 2009) suggest that the kinetic 
energy of the gas motions is at the level of 5-15% of its ther- 
mal energy in the relaxed clusters. Slight deviations from 
incompressibility should not dramatically alter our results. 

3) As shown in § 3.3, the evolution of the decay rate of 
the small-scale (/ < A„) temperature variations can be linked 
to the amount of stretching of the field lines as oc (l/s^) . Es- 
sentially the decay rate goes down because the field lines. 
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along which the heat is transported, are stretched^ The 
amount of stretching is of course limited by saturation of 
the magnetic field. We note that saturation, which can be 
one of the key effects in the problem, is completely neglected 
in our paper and is a matter of further study. Another po- 
tentially important effect we have blatantly disregarded is 
reconnection of the field lines that constantly modifies their 
topology. While we believe the simple model considered in 
this paper correctly captures the qualitative picture, direct 
numerical simulations are required to confirm this. 

4) We stress again that we only consider the suppres- 
sion of local thermal conductivity on scales I < \u- We have 
established that the fluctuating temperature gradients are 
predominantly oriented perpendicular to the magnetic fleld 
lines, which slows down thermal conduction. This process 
resembles the one that occurs macroscopically in cold fronts 
where the field lines are aligned along the cold front inter- 
face by the plasma fiow. If one is interested in the global 
heat transport on scales / 3> A„, other effects start to be 
important. In particular, the exponential divergence of the 
field lines and transverse diffusion of electrons increase the 
global conductivity (Rechester & Rosenbluth 1978; Chan- 
dran & Cowley 1998; Narayan & Medvedev 2001; Chandran 
& Maron 2004). 



5 CONCLUSIONS 

We have studied the correlations between the local fiuctuat- 
ing temperature gradients and the orientation of the frozen- 
in magnetic field lines in the turbulent ICM. One of the fre- 
quently made assumptions is that the magnetic field is ran- 
domly oriented with respect to the temperature gradients. 
We argue, instead, that strong correlation between the di- 
rection of the field and the temperature gradients can be ex- 
pected in a turbulent ICM. Gas motions tend to increase the 
temperature gradients and, at the same time, aline the mag- 
netic field lines perpendicular to the gradients. Cold fronts 
in clusters provide a vivid example of this process at large 
scales. The net result of the correlated evolution of the tem- 
perature distribution and the magnetic field is the effective 
suppression of the local heat fiux. We have calculated ex- 
plicitly the joint distribution function of the gradients and 
the angles they make relative to the field lines and demon- 
strated that significant suppression takes place for generic 
3D isotropic incompressible motions. The main results of 
this study can be summarized as follows: 

• Strong correlation of the fluctuating temperature gra- 
dients and the local magnetic fleld orientation is established 
on the eddy turnaround time. 

• For disturbed clusters, where large-scale clumps of gas 
are displaced, the largest observed gradients should be asso- 
ciated with the largest heat flux suppression. The estimates 
of the effective conductivity based on these gradients may 
not be characteristic of the bulk of the gas. 



^ The effect of the field-line stretching on the suppression of ther- 
mal conduction has previously been studied by Rosner & Tucker 
1989 and Tao 1995, but in the case of \b < ^mfp a-nd constant 
macroscopic thermal gradient. 



• On average, the decay time of temperature fluctuations 
is anti-correlated with the ampliflcation of the magnetic field 
by the gas motions. Volume averaged decay rate decreases 
with the growth of the magnetic energy density as {B'^)~^'^ . 
In other words, mixing of the ICM by the isotropic incom- 
pressible motions does not promote heat exchange as long 
as the magnetic field remains frozen. 
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APPENDIX A: STATISTICAL CALCULATION 
OF THE JOINT PDF OF fi, G AND B 

The general form of the joint PDF of the magnetic field and 
the temperature gradient is 

P(g,b,G,B;t) = {P), 

P = S(g- g{t))S(b - b{t))S(G - G{t))S{B - B{t)), (Al) 

where g, b, G and B are variables and g{t), b{t), G{t) and 
B{t) are stochastic processes that are solutions of equa- 
tions (30). Taking time derivative of P and using equa- 
tions (30), we obtain 



dtP = Lfcr^jP, 
where 

r) 

Li = -Q-ry^m ~ 9 g )9 



(A2) 






b''b')b"' 



d ■ d 

The average of equation (A2) is 

dtP = iTi^^l^P) 



(A3) 



(A4) 



and we now apply apply the Furutsu-Novikov formula (Fu- 
rutsu 1963; Novikov 1965) to calculate the right-hand side: 



{aUt)P{t)} = / dt'iaUthiit')) 






— K-/ rr. 



5P{t) 



(A5) 



where we have used equation (15). From equation (A2), 

sPit) _ I ,,, 



Sai{t) 



Lr5;c<5(t - t')p{t') 



I fmi ,,l\ dP(t ) 
+ ^i ^m[t )- 



Sai(t) 



t]P{t). 



(A6) 



The second term inside the integral vanishes by causality 
(t' < t). Using equation (A6) in equation (A5) and substi- 
tuting into equation (A4) , we arrive at a closed equation for 
the desired PDF: 



(A7) 



a p_ '^ rpij T m-r n p 
Ot-T — —Imn^i J-'j -■ • 

Since the medium is isotropic, the PDF only depends 
on G, B and the angle between the unit vectors g and b. 
Therefore, it can be factorized as 



Pig, b, G, B; t) = -^S(g^ - 1)5(6^ - l)p(^, G, B; t), (A8) 



where /i = b ■ g. The factor I/Stt^ has been introduced in 
order to keep p(/i, G, B; t) normalized to unity. Substituting 
this expression into equation (A7), we get 

IT^P = 5(3^ - l)5{b' - 1) 

x{(6'6^6™fo" -f g'g'g'^g" - g'b'g"'b'' 

-b''g^b"^g")fid^fid^ 

+{b g b g - g g'g g jfid^doG 

+ {gVg"'b" - 2b%'b"'b" + b'g'b"'g")fid^dBB 

-{g'b'g'^b" + b'g'b"'g")dGGdBB 

+g'g^g"'g"dGGdGG 

+ [2{d + l){b'b'b"'b" + g'g'g^'g") 

—2d{g b-' g b -^ b g-^b g ) — b b Oj —b-'b o^ 

-g g (<m~ g g K Jm^m 

+ [-2(d + l)/5^5'"5" + d{gVg"^b" + Vg'b^g^) 

+5V'54+5VC]aGG 

+ [-2{d + l)b'b'b"'b" + dg'b'g'^b" + db'g'b^'g" 

+b"'b"5} + Vb'^S^dBB 

+d[(d + 2){b%'b"'b" + g'g'g'^g") - d(g'y g"'b'' 

+b^g'b^g-) 

-{b'-b-S} + Vb'-Sl + g^g-Si, + /5'"C)1}P, (A9) 

where d is the number of spatial dimensions. The PDF is 
factorized, as it ought to be, and we only need to solve the 
equation for p{jj,,G, B;t). Substituting equation (A9) into 
equation (A7), we perform the convolutions involving T^„ 
[see equation (16)] using the identities 

a + 1 



T^ii^g'b'g-'b" = 



d+1 ' 



T'^ b^h^S' = 

rpij im,nsri _ (" ^ J-JV" + ^j 
^ mn^ Oj 7,-1 ' 

d+1 

2 _ 1 



T:^„b^g'b"^g" 



ning g g g 



d+1 
d-l 
d+1' 



T^ug'g"Si^ = 
T^ng'g^s^ = 



zj ^^^JX^^ ^ )\ ) 



d+1 
The result is equation (31). 



(AlO) 
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